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Abstract 

Using the Hartree-Fock-Bogoliubov factorization of the density matrix and the 
Born-Oppenheimer approximation we show that the motion of the condensate sat- 
isfies a nonlinear Schrodinger equation in the zero temperature limit. The Galilean 
invariance of the equation is explicitly manifested. From this equation some general 
properties of a superconductor, such as Josephson effects, the Magnus force, and the 
Bogoliubov-Anderson mode can be obtained readily. 

PACS#s: 74.20.-z; 67.50.Fi; 
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A large class of superconducting phenomena is well described by a static equation near 
superconducting transition temperature, proposed by Ginzburg and Landau (GL)[|TJ before 
the microscopic theory of Bardeen, Cooper, and Schrieffer||. Subsequently the GL equation 
has been derived microscopically by Gorkov|3| and extended to a wider parameter region by 
many authors^, [5]. Despite the tremendous success of the microscopic theory, the GL equa- 
tion occupies an important position, as also illustrated in the recent study of the clarifying 
of symmetry of the Cooper paring in high Tc superconductors 0, of the flux phases J/J, and 
of the insulator-superconductor transitions 0. Encouraged by the wide region of validity of 
the GL equation, several attempts to derive the equation of motion for the condensate in a 
superconductor, the so-called the time-dependent GL equation, have been made during the 
past three decadesj|, [S], [| [TIJ. Although there is no consensus for the exact form, a con- 
cern about the apparent non-Galilean invariance in the previous derived equations of motion 
exists as noted in Refs.Q, [T(J. While it can be justified near the transition temperature be- 
cause of the strong relaxation, it is questionable in the low temperature limit. On the other 
hand, the Galilean invariant form of the equation of motion has been repeatedly argued 



phenomenologically: the proposal of Ref.[TT[] in the derivation of Josephson relations; the 
proposal of Ref. || in the study of insulator-superconductor transitions, and the proposals 
based on symmetry considerations [O, O. Besides the conceptual question of the Galilean 
invariance, the different forms of the equation of motion have different physical consequences, 



as shown, for example, in the study of the decay of supercurrent |fL4[| and of the Hall effect [jTSfl . 
The requirement of the Galilean invariance will put a restriction on the form of the equation 
of motion. Furthermore, it guarantees certain global properties of a superconductor such as 
Josephson effects and the existence of the Magnus force. 

The purpose of the present paper is to give a formal microscopic derivation of the equation 
of motion for the condensate near zero temperature, which explicitly manifests Galilean 
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invariance. The present derivation is based on the existence of the off-diagonal long range 
order(ODLRO) in a superconductor [16 ] which allows a macroscopic wavefunction description 



of the condensate, and on the Born-Oppenheimer approximation which allows the separation 
of the slow motion of the condensate and the fast motion of the formation of the Cooper 
pairs. A nonlinear Schrodinger equation for the condensate is obtained. The usefulness 
of this equation is illustrated by three examples: Josephson effects, the Magnus force, and 
the Bogoliubov-Anderson mode. For definiteness, a clean s-wave superconductor will be 
considered. 

To demonstrate the generality of the present derivation and to give a firm microscopic 
ground, we first derive a set of self-consistent equations for the one-body reduced density 
matrix and the wavefunction describing ODLRO. The model Hamiltonian for N electrons 
may be taken as 

N _ h 2 1 N 

^ +5 !:</<*-■>>■ m 

Here m is the mass of an electron, r = (x,y,z), and V is the three-dimensional gradient. 
The effect of the lattice is to renormalize the electron mass and to give rise a weak, short 
range, and attractive two body interaction U with no spin dependence, which is responsible 
for superconductivity. Equation (1) is in the standard form of the Hamiltonian in the 
derivation of BCS theory. Here we ignore the effect of the electromagnetic field, and we 
shall return to it later. The motion of electrons is described by the iV-body density matrix 
D(ri 01, tn<tn', r W[, •••) t'nv'n'i t)i which satisfies the von Neumann equation 

ih^D = [H,D}. (2) 

Here a,j = ± is the spin index of jth electron. The f^-body reduced density matrix is defined 

as 

N\ N r 
ifo(ri<7i, r K a K ; r^, r' K a' K ; t) = _ J[ J d r « x 



1=K+1cji=±- 



D(rxa h r N a N ; ricr'i, r' K a' K} r K+1 a K +i, ...,r N a N ; t) . (3) 

The information needed is contained in the one- and two-body reduced density matrices. 
Then, from eqs. (2) and (3), we have the equations of motion for the one-body reduced 
density matrix R\, 

d h 2 
ift-fli( ri ±; ri±; t) = - — [V? - V?]iMn±; ri±; t) 

+ / dh 2 [U(n - r 2 ) - C/(r; - r 2 )] J R 2 (r 1 ±, r 2( x 2 ; r^i, r 2( x 2 ; t) . (4) 

Now, we assume that two-body and higher density matrices can be decomposed in 
terms of one-body reduced density matrix with ODLRO. This is the usual generalization 



of Dirac's[17| factorization which leads to the time-dependent Hartree-Fock theory. Specifi- 



cally, we have for the triplet part of the two-body reduced density matrix R 2 , 

R 2 (r 1 ±, r 2 ±; r[±, r 2 ±; t) « ^(ni; ri±; t)i?i(r 2 ±; r 2 ±; t) - i^ni; r' 2 ±; t^foi; ri±; t) . 

(5) 

Because of the existence of the existence of ODLRO, for the singlet part we have 

R 2 (n±, r 2 T] ri±, r' 2 T] t) « i2i(ri±5 ri±; t)i?i(r 2T ; r 2 ^; t) + *(n±, r 2T ; *)**(r'i±, r^; t) . 

(6) 

The existence of the wavefunction \I/ is the manifestation of ODLRO. Here, as |r — r'| — > oo, 
i?i(ri±; r^±; t) — > exponentially. In the absence of spin interaction, we have i?i(ri±; r^±; t) 
_R 1 (r 1 ; r^; t), independent of spin index, and i?i(ri±; r^; t) = 0. Therefore, from eqs. (4)-(6) 
we have for the one-body reduced density matrix Ri 

at 2m 
+ J d 3 r 2 [U(r 1 -r 2 )-U(r' 1 -r 2 ))2R 1 (r 2 ;r 2 ;t)R 1 (r l ;r , 1 ;t) 

- J d 3 r 2 [[/(n - r 2 ) - U(r[ - r^i^n; r 2 ; t)i?i(r 2 ; ri; t) 



+ J dh 2 [U{v 1 - r 2 ) - U(r[ - r 2 )]*(n, r 2 ; *)**(r'i, r 2 ; t) . (7) 
A similar factorization of the equation of motion for \1/ gives 

d h 2 

*ft-^(r 1; r 2 ; t) « -^[V? + V^(r 1; r 2 ; t) 
+ [c/(n - r 2 ) + J d 3 r 3 [f/( ri - r 3 ) + U(r 2 - r 3 )2 J R 1 (r 3 ; r 3 ; t)]] *(r 2 , r 3 ; *) 

Here \I/(ri, r 2 ; £) = [\l/(ri+, r 2 — ; t) — ^(ri— , r 2 +; t)]/\/2, and is symmetric upon interchange 
of i"i, r 2 because of the singlet paring. 

Eqs. (7) and (8) form a general self-consistent set of equations valid at any temperature. 
We note that eqs. (7) and (8) are similar to the Gorkov equations in the form given by 



Keldysh Green functions ||18|| , but here we use reduced density matrices, following the ap- 



proach of Penrose [19] in the case of superffuid helium 4. Same equations have been obtained 



previously f2(|. We emphasize that the presentation of leading to eqs. (7) and (8) is to show 
the generality of the following nonlinear Schrodinger equation for the condensate. 

For an illustration of the Galilean invariance we consider a uniform state. The uniform 
solution of eqs. (7) and (8) has the form for the one-body density matrix Ri 



fl 1 (r 1 ;ri;t) = p(r 1 -ri) 



(9) 



and for the wavefunction \1/ describing ODLRO 



*(r x , r 2 ; t) = tf (ri - r 2 ) exp [~V t\ , 



(10) 



which depends only on the relative coordinate Then it is easy to check that the uniform 
current carrying state, v = Kk/m, is also a solution of eqs. (7) and (8). In this case, 



Ri(rx; r' x ; t) = p(r x - r[) exp{ik ■ (r x - r[)} 



and 




(12) 



and eqs. (7) and (8) remain unchanged. 

To derive the equation of motion for the condensate, the theme of the present paper, we 
observe that eqs. (7) and (8) describe three kinds of distinct dynamics: the motion of the 
quasiparticle excitations, the formation of Cooper pairs, and the motion of the condensate. 
At low enough temperatures, because of the existence of the energy gap in the quasiparticle 
energy spectrum, the density of quasiparticle excitations is exponentially small. In this case, 
we only need to deal with the dynamics of the formation of Cooper pairs and the motion of 
Cooper pairs. To avoid unnecessary complications, we shall work under this restriction. 

As suggested by the demonstration from eqs. (9) and (10) to eqs. (11) and (12), we observe 
that the formation of a Cooper pair is described by the relative coordinate r = r x — r 2 . The 
motion of the condensate is described by the center of mass coordinate R = (ri + 1^2)/ 2. If 
the motion of the condensate is slow enough, which in principle can be arbitrarily slow, a 
separation of dynamics between the motion of the condensate and the formation of a Cooper 
pair is possible. The Born-Oppenheimer approximation can then be applied to get a simple 
equation of motion for the condensate, because the dynamics of the formation of a Cooper 
pair can always follow the motion of the condensate. To implement the Born-Oppenheimer 
approximation the wavefunction \l/ may be written as 



because the wavefunction for the formation of a Cooper pair is bounded in space. With 
eq. (14), the wavefunction of the condensate ip is normalized to the density of Cooper 



tt(n,r 2 ;*) = V>(R;*) 0(r,R,^,^;t) , 



(13) 



with the requirement that 




(14) 
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pairs, which is proportional to the superfTuid density. Therefore, from eq. (8) according to 



the Born-Oppenheimer approximation p2|| the wavefunction <f> satisfies following stationary 
equation 

h 2 



V{R,if),if)*) 0(r) = -— V^(r) + £/(r)0(r 

m 

+ J d 3 r 3 2i2 1 (r 3 ; r 3 ) 
+ / d 3 r 3 2i2i(r 3 ; r 3 



tf(R+E-r 3 ) + [/(R-~-r 3 ) 



[/(R+£-r 3 ) + [/(R-£-r 3 ) 



i?i(R+ -;r 3 



r ^((R-§ + r 3 )/2;*) 



^(R;t) 



(R 



r 3 J 



-i?i(R— -;r 3 



r ,V((R+| + r 3 )/2;0 



^R+2" r 3) 



(15) 



V(R;*) 

with i?x solved from eq. (7) as a function of if){if)*) and 

Given the slow varying quantities R and if)(if)*), eqs.(15) and (7) can be solved in the 
same manner as the usual Born-Oppenheimer approximation. In the special case of a uniform 



state, those equations have been solved in Ref . []2T]] (see also Ref.|20fl), where it is found that 
V = 2fi with n the Fermi energy of the system, and that \I/ determined by the gap equation 
is proportional to the electron number density n. In the more general case there will be a 
small correction to 2/x due to variation of the Cooper pair binding energy (the energy gap) 
as a function of \if)\. Now suppose that we can solve eq. (15) for V and <fi in terms of R and 
if) (if)*) generally. Using eqs. (8), (13)-(15) we find the wavefunction if) of the condensate 
satisfies the following equation: 



ih^if){R; t) = i[-iW R + a(R, if), </>*)MR; t) + [a (R, if), ^) + V{R, if), ^*)]V(R; *) 



with 



a(R, if), if)*) = J d 3 r 0* (r, R, if), if)*) V R 0(r, R, if), if)*) 



(16) 



(17) 
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and 

a (R,V,V*) = |^/ rf3r |V R 0(r,R,^V*)r-^a 2 (R,^,^). (18) 

We note that the mass of a Cooper pair in eq.(16) is twice of the electron mass, and eq.(16) 
is similar to the form in the discussion of the separation of fast and slow motions The 
the pseudo-scalar and vector potentials a and a are the response of the fast internal degrees 
of freedom, the dynamics of the formation of a Cooper pair, to the slow motion of the 
condensate. Because of the singlet pairing, <fr is spherically symmetric in a uniform state. 
In this case both a = and ao = 0. For a nonuniform state, the shape of a Cooper pair 
wavefunction is a deformed sphere, and we may choose <p to be real. In this case we have 
a = but a ^ 0. This is consistent with the generalized GL equation We point out 
that for non singlet pairing the pseudo-vector potential a may not be zero. Equation (16) 
is our desired equation of motion for the condensate, which takes the form of a nonlinear 
Schrodinger equation, and is Galilean invariant. 

The coupling to the electromagnetic field can be introduced as usual as the minimum 
coupling into the Hamiltonian eq. (1). Performing a parallel calculation, we find that the 
equation of motion of the condensate in this case satisfies the following equation: 

ih^(R; t) = -L[-ihV R - 2eA(R) + a(R, V, V*)]V(R; t) 

+ [2eA (R) + ao(R, V, *P*) + V(R, ^*)MR *) • (19) 

To reach the above equation, we have assumed in addition that the space variation of the 
scalar gauge potential A and the vector gauge potential A should be smooth over the length 
scale of the localization length of a Cooper pair to ignore the direct influence of the gauge 
potentials on the formation of Cooper pairs. This condition restricts eq. (19) to the case 
of a type II superconductor. Supplemented by Maxwell's equations for the electromagnetic 
field, the properties of a superconductor can be discussed under the above conditions. 
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Having obtained the nonlinear Schrodinger equation for the condensate, eq.(16) or (19), 
we consider in the following their consequences, Josephson effects, the Magnus force, and 
the Bogoliubov-Anderson mode. In the first two cases the dynamics of the the condensate 
is dominated by the phase of the condensate wavefunction ip. The last case is concerned 
with the compressibility of the superconducting state, involving both the the phase and the 
amplitude of the condensate wavefunction ip. 



Josephson relations. The derivation of Josephson relations in Ref. flTTf is the neatest one 



known so far, which is started by postulating of a Schrodinger equation for the motion of 



the condensate. However, no justification for the Schrodinger equation in Ref.|TT[| has been 
provided. In particular, one also needs to justify that the electric charge is 2e instead of e. 
The present derivation gives such answers. We note that by AC Josephson effect and the 
DC Josephson interference effect, the 2e of a Cooper pair is clearly established. 

The existence of the Magnus force. Since our equation of motion of the condensate, 
eq.(16), is identical to the form of Gross-Pitaevskii equation for superfluid helium 4, the 
existence of the Magnus force for a vortex line in a superconductor, which is proportional to 
the number density of the superfluid electrons, follows naturally[|15|]. This coincides with the 



result from a complete different approach |24]], therefore gives another general demonstration 



of the existence of the Magnus force. By study of the consequences, such as the Hall effect [|T5|] , 
of the Magnus force, we may have a direct measurement of the superfluid velocity therefore 
superfluid number density in comparison to the London penetration depth. 

The Bogoliubov-Anderson mode. A nonlinear Schrodinger equation may contain a phonon- 
like excitation mode|25|. We show below that this mode in eq.(16) has a propagating velocity 



of Vp/Vs with vf the Fermi velocity in the neutral superconductor, and is the plasma mode 
in the charged case of eq.(19). We write the condensate wavefunction as ip = ^/p exp{9}. 
Ignoring the small correction due to the variation of the binding energy of a Cooper pair we 
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have V = 2p(n). Expanding the phase and the Cooper pair density around their uniform 
state values, 6 = 9 + 59 and p = p + Using p oc n, and eliminating Sp in eq.(16), we 
obtain 

^^-vd^(R;t) = o. (20) 



This describes the Bogoliubov- Anderson mode, as known by Refs.g, 10 1. We note that 
the expansion leading to eq.(20) effectively chooses a reference frame. In the derivation of 
eq.(20), the relationship between the Fermi energy p and the electron number density n in 
the weak interaction limit has been used. For a real superconductor eq.(19) should be used 
together with Maxwell's equations. Because the scale gauge potential A depends on the 
superfluid electron number density, we recover the usual plasma mode. 

In conclusion, we have derived a nonlinear Schrodinger equation for the superconducting 
condensate near zero temperature, which is Galilean invariant. The Born-Oppenheimer 
approximation employed in the derivation allows a natural separation of the fast formation 
of a Cooper pair and the slow motion of the condensate. By the nonlinear Schrodinger 
equation, Josephson effects, the existence of the Magnus force, and the Bogoliubov- Anderson 
mode in a superconductor can be readily established. This shows a unified physics behind 
those phenomena. 
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